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Abstract 

The wavefronts from a point source in a solid with cubic symmetry are examined 
with particular attention paid to the contribution from the conical points of the 
slowness surface. An asymptotic solution is developed that is uniform across the 
edge of the cone in real space, the interior of which contains the plane lid wavefront 
analyzed by Burridge [6] . The uniform solution also contains the regular wavefronts 
away from the cone edge, a delta pulse on one side and its Hilbert transform on the 
other. In the direction of the cone edge the three wavefronts merge to produce a 
singularity of the form H{t)t~^^^. 



Introduction 



We consider the time dependent Green's function in solids of cubic symmetry, with focus 
on the field across the "fiat lid" associated with the conical degeneracy. Such conical 
points occur in the [111] directions in every material displaying cubic symmetry [12]. 
It is well known that the Green's function in anisotropic crystals can be split into a 
continuous field and a sequence of singularities. The latter are associated with regular 
points on the slowness surface with outward normal in the direction of the observer. 
Burridge [6] provided explicit expressions for the field contributions from regular points, 
whether they occur on sections of the slowness surface that are locally convex, concave or 
saddle shaped. He also derived the peculiar singular behavior on the fiat lids of the wave 
surface associated with the conical points. The singularity on the fiat lid is a step function 
in time, in contrast to the more singular delta pulse or its Hilbert transform (a 1/t pulse) 
which arise from regular points on the slowness surface. Furthermore, the wavefront on 
the flat lid decays with distance like 1/r^, in contrast to the 1/r decay from regular points. 
The analysis in [6] was limited to receiver positions inside the lid but not too close to 
the lid edge. Here we develop solutions that include Burridge's flat lid solution, but also 
display the transition across the lid edge. In particular a uniform solution is derived that 
incorporates the regular and the flat lid wavefronts. 
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Fig. 1. The neighborhood of a conical point on the slowness surface for iron (left) and potassium 
chloride (right). The two quasi-transverse sheets are shown for a sector of 7° angular extent 
about the Ill-axis (defined such that iC^/Csl < tan 7° in the notation of eq. (11)). 

Several studies of the elastic wavefield from conical points appeared after the original 
work of Burridge [6] . Barry and Musgrave [4] reduced the solution arising from elliptical 
conical points in tetragonal media to a single integral requiring numerical quadrature. 
Payton [13] considered the effects of double points in the slowness surface of a special type 
of transversely isotropic solid with the restriction cn = C33 = C23. This exhibits conical 
points on the symmetry axis as well as a ring of double points in the plane perpendicular 
to the axis. Based on the material restriction he was able to construct explicit integral 
solutions for a point force in the form of elliptic integrals which show that the wave front 
of conical points have lids, while the circular double points do not produce lids. Borovikov 
[5] is the only work that has considered a uniform solution valid across the edge of the 
flat lid. Borovikov's solution is, however, time harmonic and not directly applicable to the 
problem considered here, which we find more easily addressed by a suitable modification 
of Burridge's original approach. 

The purpose of this paper is to examine the way in which the singularity from the conical 
point, i.e., the flat lid wavefront, coalesces with the nearby regular wavefronts at their 
common meeting points along the directions of the conical boundary in real space. First 
we review some of the necessary preliminary results from [6] and examine the topology 
of the neighborhood of a conical point. The major results follow from the use of a locally 
parabolic approximation near the conical point, which allows certain closed form solutions 
of the Green's function. In particular, we derive a uniform approximation that includes 
the regular wavefronts and the flat lid singularity, and is uniform across the edge of the 
lid in real space. The solution also generates the singular form of the Green's function in 
the direction of the lid edge. 

The outline of the paper is as follows. The general theory for the Green's function in 
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anisotropic solids is reviewed in section 2. Materials with cubic symmetry are consid- 
ered in section 3 where the conical point is defined. The local behavior of the slowness 
and wave surfaces is discussed in section 3, and the crucial slowness surface curvature 
is introduced. The Green's function is reduced to a single integral in section 4 and the 
dominant contributions to the wave field are derived. This includes the fiat lid response 
and the regular wavefronts associated with the direct waves. Section 5 provides a uniform 
description appropriate to the conical point, using a generic function which incorporates 
the fiat lid singularity and the regular wavefront singularities. 



2 General theory for anisotropic solids 



We are concerned with the fundamental solution, G(x, t) with elements Gij, which solves 

pdlGim{y^,t) - djkidjdiGkmi^.t) = p^^^ (5(x), (la) 

Gij(x,t) = 0, t<0. (lb) 

Here S(t) is the Dirac delta and 5(x) = 6{xi)6{x2)S{x3), where the subscripts (the sum- 
mation convention is understood) i,j, etc. take on the values 1,2 and 3 and refer to a 
rectangular coordinate system. The elastic moduli satisfy the symmetries Cijki = Cjiki, 
Cijki = Ckiij, and are positive definite, Cijkiaijaki > for all symmetric non-zero a. Al- 
though there are now several alternative ways to express G(x, t), e.g. [10], we use the 
succinct Herglotz-Petrowski formula [7] derived by Burridge [6], 



G(x,t)= / dSF{^)S'{t-^-^), (2a) 
^ s 



IS 



p«)._J,^|^, ,,s. (2b) 

Here, S is the entire slowness surface; $, is the slowness vector on S\ A(^) is the associated 
unit polarization vector; the phase speed function v{$,) is a homogeneous function of degree 
one in ^, and they are related by 

Q^kik)Ak{i) = pv\l)Am. (3) 



where 

Q^k{i)=C,,kl^JL and l = i/i. (4) 

The eigenvalue equation (3) admits of three roots for w^, corresponding to the three 
sheets of the slowness surface. The integral (2a) is therefore over the union of the three 
sheets. Note that v{$,) = 1 on S", implying that the physical (dimensional) phase speed is 
f (^) = 1/^. Another property that follows from the functional form of f as a homogeneous 
function of its argument is that the gradient appearing in eq. (2b) is the wave (or group) 
velocity: c = Vt'(£). 
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3 Conical points of the cubic slowness surface 



3. 1 Elastic moduli for cubic symmetry 



We use a slightly different notation than that of Burridge [6], with the purpose of obtaining 
a higher order approximation to the slowness surface in the vicinity of the conical point. 
The elastic stiffness is 

C = 3AJ + 2//iI + 2(^2-yUi)D, (5) 

where the tensors are as follows: IS = S, JS = |(trS)I for all S G Sym, and D is a 
fourth order tensor defined by the cube axes {a, b, c}, as 



a®a®a®a+b®b®b®b+c®c 



The stiffness is positive definite if and only if /ii > 0, /i2 > and A + |/i2 > 0. The non-zero 



elastic moduli in the Voigt notation are cn = C22 
and C44 = C55 = C66 = 



C33 



c (X) c . 

2 
3^ 

A + 2/i2, C12 



(6) 



C23 



C31 = A 



The eight conical points on S occur symmetrically in each octant, although we restrict 
our attention to the conical point in the (+ + +) octant. Two illustrations of sections of S 
through the conical point can be found in Figures 2 and 3. In order to examine the local 
behavior of S we introduce the new basis set of vectors 



e.^^ (-1,-1, 2) 



e2 = -^(l, -1, 0), 



63 



73 



(7) 



When referred to these axes the moduli become, using a slight modification of eq. (3.46) 
in [3], 





C'l2 


c'l3 





^46 





'^12 




'^13 





•^46 





^13 
























C44 





^46 




^46 








C44 














'^46 





^66 



where 



c'li = A + 2/i2 



-12 



C44 



-46 



-(/i2 



(/i2 - /il). 



- 2(/i2 - /il; 



/il). 



-33 



-13 



-66 



A + 2/i2 - -(/i2 - /il) 



A + -(/i2 



/i2 - 2(/i2 - /il) 
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The acoustical tensor Q(^) defined in eq. (4) becomes in the new frame 



Q(|) =(A + 4)^ ® I + - 77(/^2 - /il) X 



2^266 ei + il + 2v^eie3 2^266 - 266 



_V2(e| - Q) - 266 2v^66 - 266 

where ^ is also referred to the basis vectors of eq. (7). 



26^ 



(9) 
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Fig. 2. The section si = S2 of the slowness sheets for Iron, where s is slowness. The inset depicts 
the conical section of half angle ^ — a. 



3.2 The slowness surface near the conical point 



Although formally exact solutions of the Christoffel equation det(Q(^)— pi) = are known 
for cubic crystals [8], we focus here on an approximation that retains the local geometry 



near the conical point. The conical point = 6^3 is a double root with 6 = y p/c'a, or 



^0 



/ 3p 

/il + 2/i2 



(10) 
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The slowness surface near the conical point may be approximated locally as an hourglass- 
type cone, such that 

C = ^o + Ce3 + 7(C,0)e0, (11) 
where the unit vector e,^ is defined 

= cos0ei + sin0e2. (12) 

The coordinates ( and < < 27r will be used to approximate the slowness surface in 
the vicinity of the conical point C = 0. 

The approximation used by Burridge [6] was 7 = f3(, corresponding to a right circular 
cone, where f3 = cot a, and ^ — a is the semi-vertical angle of the right circular tangent 
cone at the conical point, 

tan« = 4^ = v^ '^^~^^' . (13) 

0^4 yUi + 2/i2 

Note that < a < arctan\/2 = 55°, the lower limit being the isotropic value, and the 
upper limit occurs as /i2//^i 0. The angle | — a is z/ in [6]. 




Fig. 3. The conical lid is defined by the cone of semi-angle a in real space. The coordinate 
z is length along the [111] direction, and r is the perpendicular distance from the axis. The 
parameter rj is the nondimensional distance from the lid edge, eq. (30). 

The section of the cone in the plane containing the [110] and [001] axes is shown in Figure 
2. The cone in the slowness sheet is defined by the intersection of the quasi-transverse 
sheets at a point along the [111] direction. The associated wave velocity cone in real space 
is defined by the normals to the cone, which sweep out a cone of semi-interior angle a in 
real space. The fiat "lid" of this cone is the circular plane lid under investigation. We are 
particularly interested in the form of the wavefront on the lid and across its boundary 
edge (?7 = in Figure 3). The interior of the lid defines a region in space of solid angle 
(1 — cos«)27r. The eight lids of the slowness surface together define a fraction of total 
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solid angle equal to 8(1 — cos a)27r/ (47r) = 4(1 — cos a), which entails as much as 73.4% of 
space. In other words, the region of influence of the conical point resulting from a point 
source has influence in, for instance, 37.4% of space in Iron. It is important to realize that 
the field originates from a single point on the slowness surface, and therefore does not 
generate energy at infinity. The conical lid wavefront is a nearfield effect. 

In order to capture the local curvature of the slowness surface we assume a parabolic 
approximation to the curvy cone: 

7(C,0) = /5C + ^«:(0)C'. (14) 

Substituting from (14) in det(Q(^) — pi) = and expanding in powers of (, the O(C^) 
term implies (13), while the 0(C'^) term yields, in agreement with eq. (28) of [15], 



3(/ii + 2/i2 



3yUl/i2 



1 + 



'2/i2 

Kf^2 - 



9(A + /ii 



(sgn(/i2 -/ii))cos30) 



(15) 



Thus k(0) displays three-fold symmetry, with extremal values attained in the directions 
= and (/) = 7r/3. The cone angle and curvature for a range of materials are tabulated in 
Table 1. Note that the curvature for some materials changes signs, indicating that convex 
and concave sections may coexist. If /i2 — Pi > define kq = k(0), k\ = k{tt/3), and 
if yU2 ~ /^i < 0, reverse the definitions. Thus, kq > 0, and the curvature changes sign iff 
the ratio ki/kq is negative. The vanishing of k defines a parabolic line, or a line of zero 
Gaussian curvature. The examples in Figure 1 illustrate the two possibilities: the curvature 
is always positive for iron but takes both signs for KCl. Parabolic lines have been studied 
extensively, e.g. [9,15]. Shuvalov and Every [16] provide a detailed study of the shape of 
the acoustic slowness surface of anisotropic solids near points of conical degeneracy, and 
show that there can be up to three pairs of lines of zero Gaussian curvature. It may be 
checked that the condition /t = is precisely Every's condition F, eq. (16) of [9]. The 
polarization vector fields near conical points display interesting properties [6,1,2] and we 
refer the reader to the exhaustive analysis of Shuvalov [14] for details. 



3.3 The wave surface 



Let the receiver position be defined in the new coordinate system as 

X = ree + ze^, (16) 

where is orthogonal to 63 and defined by eq. (12) with replaced by 6. The contribution 
to the fundamental solution (2a) from the neighborhood of the conical point follows by 
approximating the function F of (2b) locally [6]. As mentioned above, the gradient Vf (^) 
equals the wave velocity c, and therefore |Vt'(£)| = c, the magnitude. The wave velocity 
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satisfies ^-0 = 1 [12], while near the conical point the approximations ^ ~ Co and 
^ ■ c cos a hold. The first is clear, while the second is a consequence of the fact that 
the slowness surface is to leading order conical. Hence, c ~ 1/(^0 cos a). The polarization 
has the form A = (eq. (7.15) of [6]) which together with the approximation for 

\Vv{^)\ leads to 



G 



Co cos a 
16 7r2 



dS 



Sc 



1 — cos sin 
sin 1 + cos 




5'(t-^-x) 



(17) 



Here Sc denotes the slowness surface in the vicinity of the conical point. Also, the matrix 
elements are referred to the local basis {ei, e2, e^}. It is important to note that Sc is single- 
valued and that the two distinct quasi-transverse slownesses are included by virtue of the 
fact that both C < and C > must be considered in the integral. This implies that the 
neighboring sheets at coincident values of\(\ correspond to and e^^^^ and the associated 
polarizations are (as above) A = and A = e^/2+n/2, respectively. This formulation 
obviates the need to distinguish the two almost degenerate sheets. We note that the 

surface measure is d S* = 

7^ ^ 97/90 = 



1/2 



(1 + 7^)7^ + 7^ dC d0, where 7^ = d^/d( = /3 + k(0)C and 



The difficulty in evaluating (17) arises from the phase term 

^■x = ^oZ + Cz + 7(C, 0) r cos(0 - 9), (18) 

which follows from (16) and the local approximation (11). As Table 1 demonstrates, the 
slowness surface curvature parameter k, defined in (14) varies with 0, and can change sign. 



3. 3. 1 Local approximation of the wave surface 

The wave (or group) velocity vector c is normal to the slowness surface and satisfies 
c ■ ^ = 1 . These properties, combined with the slowness in eq. (11) imply 

c = (tc^s - + (7<^/7)e^+|)/[(Co + C)7c - 7]- (19) 

In order to derive eq. (19), first note that the partial derivatives of the slowness follow 
from eqs. (11) and (12) as = 63 + 7^0^ and = 7<^e<^ -|- 7e^+|. These are tangent 
vectors to the slowness surface, and the direction of c is therefore defined by their cross 
product. This may be readily evaluated using the orthonormal properties of the triad 
{e^, e<^+|, 63}. Then using eq. c ■ ^ = 1 to fix the magnitude, we arrive at the expression 
for c. 

Consider a field point in the direction of the wave velocity: 

X = 2:63 -I- re_L = tc, (20) 
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where is the unit vector in the direction of the projection of c on the plane perpen- 
dicular to the 63 direction (the [111] axis). It follows from eq. (19) as e_L = — e<^ + 

(70/7)e0+|)/y'l + {l<j>/lY- Equating components in the directions 63 and e_L yields two 
relations, which we write as 

t = i,z + Cz- , , -= , , (21) 



These identities are exact and provide implicit relations between the field point and the 
location of the point on the slowness surface, characterized by the variables C and in 
the vicinity of the conical point. In particular, we have not yet invoked the parabolic 
approximation (14). Use of this allows us to to solve for the position on the slowness 
surface associated with a field position lying near the circular boundary of the lid in 
physical space. The edge of the lid is defined by 

2; — r/? = 0, where /? = cot a. (22) 

Table 1. The cone angle a and the nondimensionalized curvature k(0) = i<-{4>)^q/I3 follow from 
equations (13) and (15) and (3 = cot a. Here kq = k(0), ki = «;(7r/3) if /i2 — /^i > with the 
definitions reversed if — f^i < 0. Cubic elastic constants from Musgrave [12]. 



Material 


a (deg) 


Kq 


Ki/kq 


diamond 


5.32 


122.03 


1.33 


aluminum 


5.41 


117.93 


1.36 


silicon 


12.61 


22.14 


1.70 


germanium 


14.39 


16.94 


1.77 


silver chloride 


14.82 


10.61 


-0.69 


potassium chloride 


20.46 


4.29 


-3.58 


iron 


25.02 


5.11 


2.35 


nickel 


25.38 


4.93 


2.37 


gold 


28.56 


3.68 


2.87 


silver 


29.59 


3.35 


2.79 


copper 


30.95 


2.96 


2.83 


potassium 


42.16 


1.00 


4.23 


sodium 


44.84 


0.72 


4.34 


/3-brass 


45.28 


0.68 


5.10 


lithium 


53.01 


0.10 


20.33 
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The idea is to solve eq. (21)2 for ( using (14) with the assumption that z — r(] ^ 0. This 
amounts to finding solutions to leading order in the small parameter (. Thus, we find 



c 



z — rf3 



+ 0({z-rPf). 



(23) 



Substituting into (21) and again expanding in ^ gives 



{z - r(3f 



(24) 



2Kr 



The critical parameter in the previous two results is the curvature k(0), which is evaluated 
at the value of defined implicitly by eq. (20). The expansion in terms of {z — rf3) fails if 
K vanishes, and clearly a different scaling is necessary near values of (p at which fi;(0) = 0. 
We explicitly avoid such directions. 

For the remainder of the paper we adopt the major simplification that the curvature n is 
independent of and hence 7 is azimuthally symmetric. 



The removal of the variation in curvature with azimuth leads to tractable expressions for 
the Green's function, ones that are more detailed than those of Burridge [6], yet which 
retains the essential physics. At the same time, the approximation is not overly excessive 
in the sense that we expect all directions of observation will display the same general 
properties found here. The exceptions are for receiver positions defined by the possible 
parabolic lines on the slowness surface along which K{(f)) vanishes, which is outside the 
scope of the present analysis. 



4 Evaluation of the Green's function 

4.1 Reduction to a single integral 

As discussed we restrict attention to the case of k independent of 0. Under this assumption 
the azimuthal integral in (17) may be effected. First change from (p to = (p — 6, and 
noting that dS = |7(C)| a/1 + 7^ d( dip, we obtain 



7 = 7(C)- 



(25) 




1 — cos 6 cos tp sin 6 cos ip 




sin 6 cos ip 1 + cos 6 cos tp 



5'{T - R cosip). 



(26) 















10 



Here, 

T = t-^oZ-Cz, i? = 7(C)r. 
The integral can be performed using the identity 



27r 
2^ 



f{cosip) 5'{T — R cosip) 



1 d 



(27) 



(2^ 



x=T/R 



where H is the Heaviside step function {H{t) = 0, ^, 1 for t < 0, = 0, > 0, respectively) 
with the result, 



G 



^0 cos a 
87r2 



dC|7(C)l/r+7f 



T-Rcos9 RsinO 
i? sin ^ T + cos ^ 




H{R^ - T2) 

(i?2_r2)3/2- 



(29) 



Define two new dimensionless parameters 



2Kr 



T 



[z + rf3)- 



V 



r(3 



z + rf3 



(30) 



r is a reduced time variable relative to the time of arrival of the flat lid wavefront, and rj 
defines distance from the lid edge, with 77 > (< 0) if the observer lies inside (outside) 
the lid, see Figure 3. Let x = Kr(/ {z + r/3), then it may be checked that 



[z + r(3Y 
(2Kr)2 



(x+l)2-(l+r) {x - 7]f - {rf - t) 



(31) 



There are four roots io R^ — T"^ = aX x = rj ± [jf — r)^/2 qj^^ x = — 1 ± (1 + r)^/2^ 
Points near the edge of the lid are characterized by \ri\ -C 1, and the associated wavefronts 
from the lid and the "regular" wave arrive at times such that the time parameter r is 
also small. Of the four roots, all are at values of x that are asymptotically small except 
the single spurious root at x = — 1 — (1 + r)^/2 ~ _2. The origin of this root lies with the 
parabolic approximation to the sheets of the slowness surface, and it is not a real effect 
in elasticity. We therefore ignore this root and replace (31) by the cubic approximation 



R^ -T' 



,z + r(3) ^^^^^ ^.^j^ ^(a;) = (2x - r) [(x - 77)^ - (?7^ - r) 



{2Kry 



(32) 



We make the further approximation yl + (7c) " 
point. Then noting that 



CSC a, which is valid near the conical 



CT 



(i?2 - T2)3/2 



-a 



(i?2 _T2)l/2^ 
1 



(i?2_r2)3/2 r^iR^-T^/^ 



p (i?2 - r2)i/: 



2 ' 



(33a) 
(33b) 
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where 9^ is ^, etc., we deduce from eqs. (29) through (32) that 



G 



where / is the integral 



1 vr 



(3dz — cos 9 dr sin 6 dr 
sin 9 dr P dz + cos 6* 9^ 




z + r[3 



Tt ^ r ^ /X H{g{x)) 
I[T,7]) = — / dx sgn(xj — . 



(34) 



(35) 



4-2 Evaluation of the integral 



The integral (35) can always be reduced to a sum of elliptic integrals, which in turn 
depend upon the number of roots of g{x) = 0. Because g{x) is cubic it can have either 
one or three real roots. Two distinct cases must be distinguished: 



r < ?7^, three roots 
r > ?7^, one root. 

It is useful to introduce two new parameters, 

A r = r — r/^, 

which is the time relative to the curved "regular" wavefront, and 



tan$ 



Arl 



-<$<-. 

2 2 



We examine the two cases (36a) and (36b) separately. 



(36a) 
(36b) 



(37) 



(3^ 



4-2.1 Before the regular front 

We consider r < rj"^, for which there are three roots to g{x) = 0. Denote the roots by 
a > b > c, then it follows from (32) that abc = r^/2 > 0. The roots are therefore either 
all positive, or two are negative and one positive. In either case, (35) becomes 



V2 
vr 



sgn(cj 



dx 



+ 



dx 



{x — a){x — h) (x — c) {x — a){x — b) (x — c) 



(39) 
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Equations (3.131) of [11] imply that the two integrals in (39) are identical in value and 
hence they cancel if c < 0. Simplifying the integrals gives 



(40) 



TT ^Ja — c \\l a — c 
and K is the complete elliptic integral of the first kind, 

"'o V 1 - A;2 gin2 y. 

Referring to (32) we see that c > iff both r and j] are positive. Hence, replacing a, h 
and c by the appropriate roots r\ + ^Jrf — r, r] — ^rf- — r and |, respectively, yields 



(41) 



Ti^ri 2 2 
where Ar is defined in (37) and $ in (38). 



,/i (ytan($-^)), (42) 



-^^.^^ After the regular front 

We now consider the situation when (36b) holds. In this case g{x) =0 has only one root 
at X = r/2 and it follows, using equation (3.138.7) of [11], that 



TT 



dx 



1 2 

;a;-f)[(a;-r/)2 + (r-r/2)] vr V P 



where 



Substituting from (37) and (38) reduces (43) to 



IP + V-2 



2V2 



TT 



An 



'At rj^ 



-1/4 



2]) 



cosf 



" 4' 



(43) 



(44) 



(45) 



4-2.3 Summary of the two cases 

Combining (42) and (45) and using (38) provides us with a concise formula for I{T,ri): 



^ 2V2 Vcos<l> 

I = 7-; TTTT X < 

TT ArV4 



2v^sec($ - f ) K (y^tan(<l> - f )) ff(r) i/(r/), Ar < 0, 



(46) 



(cos(f - f ) 



Ar > 0. 



We are now in a position to apply these formulae and summarize the behavior of the 
Green's function. 
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5 Discussion of the wavefront singularities 



The contribution to the Green's function from the neighborhood of the conical point is 
evidently contained in the functional form of eq. (46). In addition to the flat lid response 
this should include the regular wavefronts associated with the curved surfaces of the two 
sheets near their point of intersection. 



5. 1 The flat lid wavefront 



The flrst thing we note is the step discontinuity in J at r = for 77 > 0, i.e., inside the 
lid, see Figure 3. Using (46) for Ar < and noting that $ = 7r/4 in this limit, implies 



J(0+,r/)-J(0-,r/) = — . 



(47) 



Combined with (30) and (34) this gives the contribution to the Green's function from the 
flat lid. 



G 



H(t-^oz) 



4:71 (7.2 



2 /?2^3/2 



-z —r(3 cosd r(3 smO 
r(3 sin 9 —z + r(3 cos 9 




?7 > 0. 



(4J 



This agrees with equation (7.22) of [6]. 



5.2 The curved wavefronts near the lid edge 



We next consider the singularities at Ar = for 77 < and ?7 > 0, but such that |?7| ^ 1. 
The latter implies that the fleld point lies close to the direction z = r cot a, the edge of 
the hd. 



For small r/, we have 
and therefore from (34), 



dz^ Y^dr^, - ^ Sr,, (49) 



where the unit vector A is now e_ 



G--^^AA-|^, (50) 
16 IT r-^ orj 

-e/2, 



9 9 

A = cos -ei - sin -e2. (51) 
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The behavior of / as Ar 0± depends upon how $ of eq. (38) behaves in this hmit. We 
note that 

hm $ = — sern ri. (52) 
At^o± 2 

The two cases of the observer lying just outside or just inside the hd are considered in 
sequence. 



5.2.1 The convex wavefront outside the lid 

First, if ?7 < 0, it follows from (45), (46) and (52) that 




2 . dl 



^|iJ(Ar), -^2^\6ir-v'). (53) 

The behavior of the Green's function near to but outside of the lid edge is therefore, using 
(30), (50) and (53), 



G.M!AA-ViN4<-J..-(i^). ,<0. (54) 

47r Kr \ Zur J 

The delta function singularity is what one would expect from the locally convex shape 
of the slowness surface and the fact that k > is assumed, as occurs, for instance, for 
iron in Figure 1. Equation (54) also agrees with the general result of eq. (6.8) of [6]. It is 
interesting to note that the strength of the singularity decreases as as the field point 
approaches the lid edge. 



5.2.2 The saddle-shaped wavefront inside the lid 

We next consider the singularity at Ar = for 77 > 0, for which we use 



1 



K{q)^\n^=^, q^l. (55) 



Equations (42) and (52) together imply that for fixed but asymptotically small and posi- 
tive rj and Ar = 0(77), 

2 1,1. 
/^--=ln-— , Ar^O±. (56) 



TT \^2rj I Ar| ' 
Hence, 

dl 2 ^/2Tj 



dl] 71 At + iO ' 

and the Green's function is 



(57) 



G'^l^AA^^ \.^,y.. , ,>0. (68) 

47r Kr nU- ^^z - ^ + ^0 1 
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Again, this has the expected form for a singularity associated with a part of the slowness 
surface which is saddle shaped, in agreement with equation (6.9) of [6]. We note that the 
magnitude of the singularity also depends upon y^. 

The difference between this case and the prior result eq. (54) is that the the normal to the 
slowness surface in the latter case is exterior to the cone, while for the present situation 
it lies inside the conical surface. With reference to the picture for iron in Figure 1, the 
convex wavefront arises from the lower convex part, while eq. (58) is associated with the 
interior of the conical "bowl". 



5.2.3 Summary of the response inside and outside the lid 

The three singularities associated with the lid wavefront and the wavefronts from the 
convex and saddle-shaped regions near the conical point are described by equations (48), 
(54) and (58), respectively. In summary, from (54) and (58): 



G 



S {t - 



{z-rl3y 
2Kr 



r]<0, 



r]> 0. 



(59) 



These results are, so far, predictable from the analysis in [6], combined with the parabolic 
form of the slowness surface assumed here. In that sense we have not yet derived any 
new forms of the Green's function, other than to specialize them to the parabolic local 
approximation of the slowness surface. The remainder of the paper generates expressions 
that are not special cases of prior results. We conclude this section by considering the 
field right on the lid edge. The transition from one regime to another is encapsulated in 
the general expressions of equations (46) and (50) and is examined next. In section 6 we 
derive uniformly asymptotic expressions that are valid across the lid edge. 



5.3 The singularity in the direction of the lid edge 



The lid edge direction z = r cot a is particularly interesting, and the singularity there 
corresponds to the origin of the nondimensional coordinates: r] = t = 0. The general 
expressions in (46) may be used. In this limit the domain of the term for At vanishes, 
implying zero field for r < 0. The field after r = follows from (46) and (50). Using (38) 
it is straightforward to show that 



drj 



H(t) 



r]=0 



7r(r + ^ 



1 dK 



dk \^/2 




(60) 



The derivative of the elliptic integral in (60) simplifies to tx j \plK for the particular argu- 
ment 1/^2 (see eqs. (8.123.2) and (8.122) in [11]). Then using the identity K{\l^/T) = 



16 



r^(l/4)/4y^ from (8.129.1) of [11], where r(l/4) ^ 3.62561, and making the replacement 
1 + r/4 — >• 1 in (60), we obtain 



di_ 

dr] 



H(t) 2V2^ 



r3/4 r2(l/4) 
The Green's function is therefore, from (50), 



G 



in 



A A 



T 



(2«;)V4 H{t-^oz) 



(3 r2(l/4) nr[r{t-^ozYY/^ 



167r2 



T\l/A) (K{t-^^z) 



(61) 



1/2- 



167r2 



2(3z 



ri = 0. 
(62) 



6 A uniform solution in rescaled variables 



The various singularities and their interaction can be understood by a rescaling suitable 
for the neighborhood of the lid edge. Let 



T = ST] 



(63) 



where the time-like parameter is 



2Kr 



{z - r/5)2 



(64) 



We recast the results in terms of s and rj, where the region of interest is near the lid edge, 
i.e., |?7| <C 1 but s = 0(1). The three wavefronts discussed in the previous section, the 
flat lid wavefront, and the "regular" wavefronts inside and outside the cone z — /3r = 
correspond to (s = 0,?7 > 0), (s = 1,?7 > 0) and (s = 1,?7 < 0), respectively. The rescaled 
time allows us to separate the three, and to observe their coalescence. 



Since |?7| <^ 1, (38) implies that tan$ ^ sgn t]/ ^\s — 1|, and after some simplification, 
(46) reduces to 

2 

/(r,?7) ^ ^= J(s,sgn?7), (65) 



where 



J{s, sgnr/j = X < 

TT 



s < 0, 



2 (l + v^r^) K (y^^) H{sgnr]), < s < 1, 



(66) 



-1/4^ 



1 _|_ sgn rj \ 



S>1. 
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Note that J only depends on rj via its sign. It follows that 

dl — sgn 7] 



dr] 



\V 



13/2 



7 . 9J 

OS 



for ?7 7^ 0, 



(67) 



and so the fundamental solution in the vicinity of the cone z = r cot a (although not on 
it, i.e. ?7 7^ 0) may be written, from (50) and (67), as 



An 



^2 1 3/2 



J{s,r]) + As^{s,r]) 



(68) 



2.5 




-1 1 2 3 4 5 

S 



Fig. 4. The function J(s, sgn ry) of eq. (66) 



The function J is plotted in Figure 4. Note the step singularity at s = for 77 > 0, 
corresponding to the lid wavefront of (48). The singularities at s = 1 are a step for 
?7 < 0, which results in the delta pulse of (54), and a logarithmic singularity for 77 > 0, 
corresponding to the 1/t wave in (58). 

The propagating singularities are all contained in the expression [J + AsdJ/ds) in (68). 
Because this can exhibit delta function behavior which is difficult to display on a plot, we 
look instead at the response from a step load, i.e., the solution to (la) with 5{t) replaced 
with H{t). Integrating (68) gives the step response near the edge of the lid as 



An 2K,r 

where 



^AA^^M(s,sgnr/), (69) 



rs / gj\ 
M(s, sgn?7) = — sgn ?7 J \J + As — jds 

= sgn 7] (3 J^' J ds - As (70) 

The function M is plotted in Figure 5. For observers outside the cone, i.e. z — rf3 < 0, 
there is no signal until the arrival of the regular wavefront at s = 1 (r = 77^), and is the 
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step function instead of the delta response in eq. (59). Inside the cone [z — r(3 > 0) the 
first, weak, signal is from the lid edge at s = (r = 0) which has the form of a ramp 
function, as compared with the step in eq. (48). This is followed by the regular wavefront 



at s 



7] ) with a logarithmic singularity. 



M 



-2 



-4 











ri < 






X"^ ri >0 



0.5 



1.5 



2.5 



Fig. 5. The function M(s,sgn7y) of eq. (70) 

The one remaining feature not evident in the uniform solution in terms of the J and M 
functions is the behavior right on the lid edge. This follows by taking the limit s ^ oo 
with rj ^ 0. Using (66) and integrating by parts, gives 



M(s, sgn 77) — M(l, sgn 77) 




P(l/4) 



(71) 



Thus M(s, sgn 77) becomes independent of sgn 77 for large s, and grows as 1.5255 which 
is consistent with the singularity along the directions for which 7/ = 0, see equation 

(62). 



7 Conclusion 



We have revisited the problem of determining the waves emanating from the vicinity of 
the conical point in materials with cubic symmetry. The analysis is based on the general 
formulation of Burridge [6] and employs a higher order approximation to the slowness 
surface. The key ingredient is the parabolic approximation (14) for the vicinity of the 
conical point, which depends upon the curvature k of eq. (15). By assuming the curvature 
is constant in the azimuth the evaluation of the Green's function is reduced to a single 
integral (35). This reproduces the flat lid singularity (48) derived by Burridge [6] as well 
as the regular wavefronts originating near the conical point. The wavefront singularity on 
the lid edge, eq. (62), has been evaluated. A uniformly vahd representation is presented 
for the transition as the observation direction crosses the lid edge. This involves a single 
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function, based on rescaled space and time coordinates, that displays the regular wavefront 
singularities and the flat lid singularity. 
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